Blistering of viscoelastic filaments 
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When a dilute polymer solution experiences capillary thinning, it forms an almost uniformly 
cylindrical thread, which we study experimentally. In the last stages of thinning, when polymers 
have become fully stretched, the filament becomes prone to instabilities, of which we describe two: 
A novel "breathing" instability, originating from the edge of the filament, and a sinusoidal instability 
in the interior, which ultimately gives rise to a "blistering" pattern of beads on the filament. We 
describe the linear instability with a spatial resolution of 80 nm in the disturbance amplitude. For 
sufficiently high polymer concentrations, the filament eventually separates out into a "solid" phase 
of entangled polymers, connected by fluid beads. A solid polymer fiber of about 100 nanometer 
thickness remains, which is essentially permanent. 

PACS numbers: 47.20.Dr, 47.20.Gv, 47.57.Ng 



When a drop falls from a faucet, surface tension drives 
the fluid motion toward breakup in finite time, and a drop 
separates. This pinch-off occurs in a localized fashion 
[l[ , and the neighborhood of the point of breakup is de- 
scribed by a similarity solution Q . If however very small 
amounts of high molecular weight polymer are added, 
an almost perfectly cylindrical thread is formed instead 
@, 0, S @| ■ The reason is that wherever there is a lo- 
cal decrease in radius, fluid elements are stretched, and 
the polymers along with it. This will increase the ex- 
tensional viscosity of the fluid-polymer mixture 0|, and 
further flow is inhibited, thus forming a uniform and sta- 
ble filament. 

For most of this paper, we produce a filament by plac- 
ing a drop of liquid between two solid plates, which are 
rapidly drawn apart Q. (In a simple and educational 
version of this experiment, a drop of saliva is placed be- 
tween thumb and index finger). A single filament forms 
between the plates, which thins as surface tension drains 
fluid from the filament, and into two roughly hemispher- 
ical reservoirs at the endplates. 

The present paper addresses the later stages of the 
thinning of the polymeric filament, when polymers have 
come close to their full extension. Thus the mechanism 
that formerly used to stabilize the thread is no longer ef- 
fective, and tiny beads begin to appear on the filament, 
(see Fig. [TJ images 5 and 6) |{|. We will refer to this 
process central to the present study as "blistering" . This 
instability occurs when the filament is only \1\xm in ra- 
dius, requiring extreme spatial resolution. If the concen- 
tration of polymer is above 1000 ppm, the filament can 
become very long-lived compared to the timescale of a 
dissolved polymer. 

Theoretically, the period of exponential thinning has 
recently been described within a long- wavelength descrip- 
tion Nonetheless, the full three-dimensional, axisym- 
metric problem remains unsolved. The effect of finite 
polymer extensibility has been studied numerically in 
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FIG. 1: Color online. The minimum radius hmin as a function 
of time. Not all data points are shown. Between number 1 
and 3, the curve is well described by h min (t) — hoexp(—t/r) 
with r = 130 ± 30ms (red, straight line). A plateau is reached 
between 3 and 4 which is associated with an instability near 
the endplates, followed by rapid pinching. Between 5 and 6 
the curve can be approximated with two linear laws (dotted 
blue and dashed green lines). For further explanation see text. 



[icj |. once more using a long- wavelength model. The fila- 
ment is found to fail near its end via a localized similar- 
ity solution, in contrast to the much more complex sce- 
nario found here. The first clear experimental description 
of blistering is found in [ll| , which focuses on the later 
stages of the instability, in the course of which droplets 
with a hierarchy of sizes are found. We first focus on the 
onset of the blistering instability, for which widely di- 
verging theoretical exp lanations have been expressed in 
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Our experiments were performed with aqueous solu- 
tions of polyethylene oxide (PEO) of molecular weights 



M w = 1 — 8x 10 6 amu and concentrations 0.010 < c < 0.2 
weight percent. While the formation of a filament and 
its subsequent instability could be observed well below 
the overlap concentration Q of polymers, we focused on 
higher concentrations, as the process is slower and easier 
to observe. Our reference system has Mw = 4 x 10 6 amu 
and c = 0.2 weight percent, with c ov = 0.07 weight 
percent the overlap concentration. The samples were 
characterized with a Thermo Haake MARS rheometer 
using cone plate geometries. The zero shear viscos- 
ity was r]o ~ 50mPas and in the range of shear rates 
0-1 < 7 < 2000 shear thinning was present down to a 
value of rjoo w 4mPas. The surface tension was deter- 
mined by the pendant drop method to 7 ~ 60.9mN/m 

H El- 

We used a capillary breakup device similar to the one 
described in [ll|, [l5|. This setup is also commercially 
available to measure the extensional rheology of suspen- 
sions (CaBER, Thermo Fisher Scientific, Karlsruhe, Ger- 
many). To ensure maximum reproducibility, we used the 
following protocol: plates of diameter d — 2mm were 
held at a distance of I = 2.5mm for the purpose of relax- 
ation for several seconds. Then the plates are drawn to 
I = 3.5mm within 40ms, only slightly exceeding the limit 
at which a capillary bridge of the solvent experiences a 
Raylcigh-Platcau instability and breakup. The thinning 
process is observed with an IDT X- Vision X3 digital high- 
speed video camera with Nikon Microscope objectives of 
up to 20 x magnification. At the highest magnification, 
the diffraction limited resolution is 0.6/xm, and depth of 
field is 5/im. A Halogen back-light allows frame rates up 
to 6000 frames per second and exposition times down to 
10/is. 

Figure [1] shows a typical recording of the thread radius 
h m in(t) in the cylindric region in a semi-logarithmic plot. 
For plug flow in a cylindrical filament, the elongation rate 
is determined from e = —2dlnh/dt, thus e is constant for 
most of the filament thinning, which follows an exponen- 
tial law (the regime between 1 and 2 in Fig. [T]) [4]. The 
axial stress a zz supported by the polymers balances the 
increasing capillary pressure j/h m i n , which means that 
the extensional viscosity tje = Czz/t = l/{h m in^) a ls° in- 
creases exponentially. At 3 (see Fig. Q]) the thread radius 
first has a plateau, and then thinning accelerates again. 
The reason for the acceleration is that the polymers have 
almost reached their maximum extension, so their exten- 
sional viscosity can no longer increase. But this means 
that e has to increase rapidly, implying a steep increase 
in the slope of \ogh(t), as seen in Fig. [TJ Once t]e has 
reached a plateau, which we estimate at h m i„ « 12/im 
to be t]e (12/im) w 330Pas, the filament behaves essen- 
tially like a Newtonian fluid [l(j, and is thus subject to a 
capillary instability [l6| . as confirmed below. Note that 
this value of the extensional viscosity corresponds to an 
increase by 5 orders of magnitude over r) wa t er = 10 Pas 
of the solvent. 
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FIG. 2: Growth of a sinusoidal instability of the viscoelastic 
filament that develops into a group of droplets on the thinning 
filament. The spacing of the pictures is 300~ 1 s. The time 
window between 4 and 5 in Fig.[Tjis represented by the images 
up to 0.287s showing the range of exponential growth. 
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FIG. 3: Main: The growth of the amplitude of a sinusoidal 
surface deformation on a filament of radius Rq — 10/im. The 
origin of the time axis has been shifted relatively to Fig. [TJ 
Circles and crosses are experimental data of two different 
runs, from the latter we were able to detect amplitudes as low 
as 80 nm. The straight and the dotted line are exponential 
fits, giving an inverse growth rate of uj — 9.3 ± 0.1ms. Inset: 
The sinusoidal surface deformation at t = 18.7ms. Points are 
experimental data and the solid line is a fit with a sine and a 
linear offset. The selected wavelength is X/Ro = 12 ± 0.9. 



Before we describe the novel instability that occurs be- 
tween 2 and 3 in Fig.Q]), which is localized near the end- 
plates, we concentrate on the subsequent spatially uni- 
form, linear instability which is shown in Fig. [5] (between 
4 and 5 in Fig. [I}. At first, no oscillations are visible on 
the images of Fig. [2 however, as seen in Fig. [3l we are 
able to resolve perturbations down to an amplitude of 
A = 80nm, corresponding to significant super-resolution 
17| . This is done by fitting the profile with a sine func- 
tion with wavenumber, phase and amplitude as free pa- 
rameters over many wavelengths (see inset). The algo- 
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FIG. 4: a) The final state of the filament. Beads are formed 
off-center relative to the thread, b) Scanning Electron Mi- 
croscopy image of two beads, connected by a thread (interme- 
diate resolution). The structure was caught and dried upon 
the substrate c) Another example of the structure, the red 
box indicating a closeup at high magnification shown in d). 
The diameter of the fiber can be as small as 70nm. 



rithm converged down to the stated maximum resolution. 
The last four pictures of Fig. [2] show the beginning of the 
non-linear stages of the instability, finally leading to the 
formation of smaller secondary droplets [ll|, [l3| ■ In the 
main panel of Figj3] we plot the growth of the sinusoidal 
approximation over time, an example of which is shown 
in the inset. Over more than a decade, the growth is 
very well described by an exponential, providing a clear 
signature of a linear instability, which develops uniformly 
in space. 

/From a fit to the exponential, we find an inverse 
growth rate of l/u> = 9.3 ± 0.1ms. Linear stability of 
a viscous fluid thread [3] predicts cj = ~//(6RoT) e ff), 
which gives an estimated extensional viscosity of rj e f f = 
9Pas±2, more than one order of magnitude smaller than 
the extensional viscosity t\e (12/Ltm) estimated above. 
At the same time, we are able to fit - as expected 
[l|| - a linear law h m in = —0.44 x 10~ 3 m/ sAt in the 
range 8^tm > h m % n > ifj,m. Comparing to the law 
hmin = 0.07"f/r) e ff At for viscous pinching [18|, this gives 
r] e ff = lOPas, which agrees nicely We do not have 
a ready explanation for the discrepancy between 7y e // 
and rjE (12/ito). Among possible explanations are non- 
uniformities in the polymer concentration (see below), 
and transient relaxation of polymer stresses during the 
plateau between 3 and 4 in Fig. [T] when there is no flow. 
A kink in the linear shrinking behavior at h m i n ~ 3.8/ito 
toward a less steep slope of —0.17 x 10 _3 m/.s could cither 
be seen as a first indication for the onset of a draining 
process discussed below, or as as a transition from the 
viscous to the inertial- viscous pinch-off regimefljj]. 



The formation of successive generations of beads has 
already been studied extensively 11 , 13| ; we focus on the 
very final stages of the thinning process, when the for- 
mation of new beads has come to rest (6 in Fig. [T|) . If the 
polymer concentration was greater than 1000 ppm, the 
filament connecting two beads never breaks, and a pat- 
tern as shown in Fig [4^) is formed. What is remarkable 
is that most beads are off-center with respect to the fila- 
ment. Comparison with the problem of fluid drops on a 
fiber [l^l shows that there must be a finite contact angle 
between the drops and the filament for such a symme- 
try breaking to occur. In other words, the thin filament 
must have formed a (solid) phase different from that of 
the drops. 

To confirm this idea, we produced the Scanning- 
Elcctron-Microscopy (SEM) images shown in Fig.|U Ob- 
ject slides were pulled quickly through the liquid bridge, 
before the accumulated elastic stress would lead it to 
retract into a single droplet. Panel b) shows two rem- 
nants of two droplets being connected by a persistent 
thin thread. Increased magnification (panels c) and 
d)) allowed us to estimate the diameter of the fiber as 
75 — 150nm. Assuming a constant polymer concentra- 
tion of the solution, and taking for the fiber the density 
of PEO, the amount of polymer in such a fiber equals 
a fluid diameter of 3/Ltm, thereby representing a lower 
bound for the onset of the concentration process that 
leads to this solid fiber, but which is likely to start ear- 
lier. Our physical picture is that polymers become entan- 
gled, while solvent drains from the filament, leading to 
even higher polymer concentration and increased entan- 
glement. Further evidence for this concentration process 
was already found in [2(|, where birefringence measure- 
ments were performed to examine molecular conforma- 
tions in the break-up process. We can exclude evapora- 
tion to be a factor in the formation of solid fibers, based 
on our estimates of evaporation rates, as well as prelimi- 
nary experiments in a two-fluid system. 

Finally we would like to describe the transition regime 
(between 2 and 3 in Fig. []} where an instability of the 
homogenous elongational flow originates from the bound- 
ary, and is reminiscent of phenomena reported in [l2T ]. 
The transition region at the edge of the filament was de- 
scribed in detail in Q for the case that polymers are far 
from stretched. At full stretch, the transition region con- 
stricts (as seen in the last shape in the inset of Fig. [5] 
(full triangles)), thus inhibiting the flow out of the fil- 
ament. The most sensitive probe for the flow are tiny 
bubbles inside the filament, whose trajectories arc shown 
in Fig. [5] During the constriction phase, the flow stops, 
and the tracer positions form a plateau. In the absence 
of flow, polymers relax inside the filament, and elastic 
stresses are eventually no longer able to sustain the cap- 
illary pressure. Fluid from the filament empties into the 
end cap, causing a sudden flow, which appears as a jump 
in bubble position in Fig. [5] The process repeats itself 
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FIG. 5: Trajectories of two bubbles in the filament; they move 
in parallel, indicating uniform flow inside the filament. First 
the bubbles are convected by the extensional flow produced 
by the thinning. Then their position undergoes a sequence 
of plateaus, associated with consecutive contractions and re- 
laxations - the "breathing" at the very end of the filament, 
shown in the inset. The breathing correspondents to the red 
shaded region between 2 and 3 in Fig. [1] At 3, the bubbles 
run out of the field of view. 

periodically, on a timescale that increases from step to 
step, but which is of the same order of magnitude as the 
polymer relaxation time. 

The height of the final jump gives a characteristic 
length scale of 150^m, which is comparable to the wave- 
lengths of periodic disturbances on the filament (FigEJ). 
In principle, each of the plateaus shown in Fig. [5] should 
result in a corresponding plateau in h m i n . However, the 
plateaus are too small to be resolved, apart from the 
last two between 3 and 4 in Fig. Q] At the end of the 
plateaus the filament shrinks again by draining liquid 
into the reservoir at the endplates or into a large bead 
that typically forms in the middle of the filament, see 0. 

In conclusion, we have demonstrated three key phe- 
nomena: 

(i) Between 2 and 3 in Fig.[TJ the trumpet-shaped transi- 
tion region connecting the filament to the reservoirs con- 
stricts periodically, interrupting the flow. 

(ii) The blistering instability of the filament, which leads 
to beads, is a linear capillary instability (between 4 and 5 
in Fig. QJ. As the polymers reach full stretch, their con- 
tribution is once more Newtonian, but with a viscosity 
that is many times that of the unstretched state. 

(iii) Using electron microscopy, we provide evidence that 
the filament remains intact in its latest stages, because 



polymer strands become sufficiently concentrated to be- 
come solid-like. The most compelling evidence for this 
fact is that the contact angle between the thread and the 
beads sitting on it becomes finite. Existing theoretical 
models are clearly inadequate in addressing this polymer 
behavior at full stretch. 
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